LMU-TPW 1996-29 and LBNL-39570 



Quantum Symmetry of Hubbard Model Unraveled 
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Superconducting quantum symmetries in extended one-band one-dimensional Hubbard models are shown to orig- 
inate from the classical (pseudo-)spin symmetry of a new class of models; the standard Hubbard model is a special 
case. The quantum symmetric models provide extra parameters but are restricted to one dimension. All models 
discussed are related by generalized Lang-Firsov transformations, some have symmetries away from half filling. 
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The exploration of high temperature superconductiv- 
ity in cuprates has greatly revived interest in the Hub- 
bard model |jl| as a model of strongly correlated electron 
systems [||-Q. Despite its formal simplicity this model 
continues to resist complete analytical or numerical un- 
derstanding. Symmetries of the Hubbard Hamiltonian 
play a major role in the reduction of the problem. They 
have for instance been used to construct eigenstates of 
the Hamiltonian with off-diagonal long range order H , to 
simplify numerical diagonalization H and to show com- 
pleteness of the solution |0] to the one-dimensional model 
[^ . Since the work of Heilmann and Lieb Q the Hubbard 
model is known to have coupling constant independent 
symmetries (spin and pseudospin) and coupling constant 
dependent "hidden" symmetries in addition to spatial 
symmetries for special configurations (benzene-hexagon 
in ly]). Operators corresponding to the hidden symme- 
tries were found in [ |lO|Jll| ; see also the nice review III4 . 

The (pseudo-)spin SO (4) symmetry [||Jl|l of the 
standard Hubbard model is restricted to the case of 
an average of one electron per site (half- filling) , so re- 
cent speculations jlj] about extended Hubbard models 
with coupling constant independent generalized (quan- 
tum group) symmetries away from half-filling attracted 
some attention. A careful analysis of the new models 
reveals that this quantum symmetry exists only on one- 
dimensional lattices and in an appropriate approximation 
seems still to be restricted to half-filling. Despite these 
shortcomings the existence of novel symmetries in Hub- 
bard models is very interesting and worth investigating. 
Quantum symmetries of the Hubbard model were first 
investigated in the form of Yangians |15[| ; quantum su- 



persymmetries have also been considered [ [l6[ . 

In this letter we shall investigate the origin of quantum 
symmetries in extended Hubbard models. We will find 
a one-to-one correspondence between Hamiltonians with 
quantum and classical symmetries. Guided by our results 
we will then be able to identify models whose symmetries 
are neither restricted to one-dimensional lattices nor to 
half filling. 

Originally introduced as a simplistic description of nar- 
row d-bands in transition metals, the Hubbard model 
combines band-like and atomic behavior. In the stan- 
dard Hubbard Hamiltonian 



^Hub = " X! 
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this is achieved by a local Coulomb term and a competing 
non-local hopping term. Here aj^, Uia are creation and 
annihilation operator^ for electrons of spin a £ {1,1} 
at site i of a D-dimensional lattice, (i,j) denotes nearest 
neighbor sites and n^o- = aj^aicr- The average number of 
electrons (^^ ^ Uiu) is fixed by the chemical potential ^. 
The standard Hubbard model has a SO (4) symmetry at 
/i = m/2, the value of /i corresponding to half filling in 
the band-like limit. This symmetry is the product of a 
magnetic SU(2)j^ (spin) with local generators 



^m — ^t^ij 



X^ — a\a 



i"T' 



Hm. ^ nt - ni 



(2) 



^We will use the convention that operators at different sites 
commute. On a bipartite lattice one can easily switch to an- 
ticommutators without changing any of our results. 



and a superconducting SU(2)^ (pseudo-spin) with local 
generators 

X+ = a|a|, X~ = a^a|, 77^ = n| + n^ - 1, (3) 

modulo a Z2, generated by the unitary transformation 
{a I 'r-^ a.) that interchanges the two sets of local gen- 
erators. The mutually orthogonal algebras generated by 
(^ and (||) are isomorphic to the algebra generated by 
the Pauli matrices and have unit elements Is = -ff^, 
Im = H^ with Is -|- Ifri = 1- The superconducting gen- 
erators commute with each term of the local part ij('°^) 
(first two terms) of the Hubbard Hamiltonian (|l|) pro- 
vided that /i = u/2. This can either be seen by direct 
computation or by studying the action of the generators 
on the four possible electron states at each site. It is also 
easily seen that the magnetic generators commute with 
each term of i/('°'^) ; in the following we will however focus 
predominantly on the superconducting symmetry. 

To check the symmetry of the non-local hopping term 
we have to consider global generators O: These genera- 
tors are here simply given by the sum ^ Oi of the local 
generators for all sites i. The rule that governs the com- 
bination of representations for more than one lattice site 
is abstractly given by the diagonal map or coproduct A of 
U{su{2)). Generators for two sites are directly obtained 
from the coproduct 



\X+,X-] = ^^ '^ " , \H,X^] =. ±2X±. (4) 

L J q — q ^ '- -' 

(The proof uses H^ = Hg.) It immediately follows that 
^(loc) j^g^g ^ local quantum symmetry. As is, this is a 
trivial statement because we did not yet consider global 
quantum symmetries. Global generators are now defined 
via the deformed coproduct of SUg(2), q G R\{0} 



A,(X±) - X± ® g-^/2 + q"/^ ® X-- 



(5) 



The local symmetry can be extended to a non-trivial 
global quantum symmetry by a modification of the Hub- 
bard Hamiltonian. The idea of [H was to achieve this 
by including phonons. Before we proceed to study the re- 
sulting extended Hubbard Hamiltonian iJextj we would 
like to make two remarks: (i) We call a Hamiltonian 
quantum symmetric if it commutes with all global gener- 
ators. This implies invariance under the quantum adjoint 
action and vice versa, (ii) Coproducts of quantum groups 
are coassociative but not cocommutative. This means 
that the reduction of global symmetry to that of next- 
neighbor terms holds only for one-dimensional lattices. 
The practical implication is an absence of quantum sym- 
metries for higher-dimensional lattices. (For a triangular 
lattice this is illustrated in Fig. ||) 



A(X±)=X±®1 + 1®X±, A{H)=H®l + l®H, * ~ \~ / ^ A is cocommutative 

• • • • • • 



while generators for N sites require {N — l)-fold iter- 
ative application of A. Coassociativity of A ensures 
that it does not matter which tensor factor is split up 
at each step. Another distinguishing property of this 
classical coproduct is its symmetry (cocommutativity) . 
This property and coassociativity ensure that we can ar- 
range that the two factors of the last coproduct coin- 
cide with any given pair {i,j) of next-neighbor sites; see 
Fig. |l|. It is hence enough to study symmetry of a single 
next-neighbor term of the Hamiltonian to prove global 
symmetry. 

(id (g) ■ ■ ■ (g) id <E) A (E) a g) ■ ■ ■ g) id) 



FIG. 2. In _D 7^ 1 symmetry of next-neighbor terms implies 
global symmetry only if A is classical. 

The extended Hubbard model of Km (with some mod- 
ifications pq ]) introduces Einstein oscillators (parame- 
ters: M, Lu) and electron-phonon couplings (local: A- 
term, non-local: Y[a,Tij^): 



-ffnxt = u 



^ n^in^i - ^ ^ n,<j - A • ^ 
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Y^ a'l„a^aT^Ja, (6) 
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with hopping amplitude 



FIG. 1. Coassociativity of A reduces global symmetry to 
symmetry of next-neighbor terms {i,j) if D = 1. 

The search for quantum group symmetries in the Hub- 
bard model is motivated by the observation that the local 
generators X^, X~ and Hg in the superconducting rep- 
resentation of SU(2) also satisfy the SU^(2) algebra as 
given in the Jimbo-Drinfel'd basis [E7| 



Tija = TJ^a = t exp(Cey • (x^ - x/) + in ■ {pi - pj)). (7) 

The displacements Xi of the ions from their rest positions 
and the corresponding momenta pi satisfy canonical com- 
mutation relations. The e^ are unit vectors from site i to 
site j. For H —Q the model reduces to the Hubbard model 
with phonons and atomic orbitals "(/"(r) ~ exp(— (^r) in s- 
wave approximation [n8[. 



The local part of iJoxt commutes with the generators 
of SU, (2)^ iff 



A2 



/i = 



Mw2 



(8) 



(For technical reasons one needs to use modified genera- 
tors X^ = exp{^ ^l,,'^2 )X^ here that however still sat- 



" umlo^ , 



isfy the SUg(2) algebra.) 

The nonlocal part of i/cxt and thereby the whole ex- 
tended Hubbard Hamiltonian commutes with the global 
generators iff 



A — hMuP'K, q — ex.p{2KCh), 



(9) 



where k 



iij ■ K for i,j ordered next neighbour sites. 
For q =/= I the symmetry is restricted to models given on 
a 1-dimensional lattice with naturally ordered sites. 

From what we have seen so far we could be let to the 
premature conclusion that the quantum symmetry is due 
to phonons and that we have found symmetry away from 
half filling because fj, ^ u/2. However: the pure Hubbard 
model with phonons has k = and hence a classical sym- 
metry {q = 1). Furthermore X ^ implies non- vanishing 
local electron-phonon coupling so that a mean field ap- 
proximation cannot be performed and we simply do not 
know how to compute the actual filling. Luckily there is 
an equivalent model that is not plagued with this prob- 
lem: A Lang-Firsov transformation with unitary opera- 
tor U = exp(iK • J2j Pj'iT'ja-)- leads to the Hamiltonian 



H„ 



UH,^tU~^ = i/ext(A',M',Ai',T',ja), (10) 



of what we shall call the quantum symmetric Hubbard 
model. It has the same form as i?cxt, but with a new set 
of parameters 



A' ^X~ Mcj^hK 

u' = u- 2h\ ■ K + Moj'^h^K^ 



^! — ^i + hX- K. 



1 



-Muj^h^K^ 



(11) 

(12) 
(13) 



and a modified hopping amplitude 



TL 



ty (1 + {q~ - l)n,a){l + [q— - l)nj,) (14) 



where iy = iexp(Ceij • {xi — Xj)). The condition for 
symmetry expressed in terms of the new parameters is 

A' = 0, /i' = ^, (15) 



u 



i.e. requires vanishing local phonon coupling and cor- 
responds to half fillingl tij may also be turned into a 
(temperature-dependent) constant via a mean field ap- 
proximation. This approximation is admissible for the 
quantum symmetric Hubbard model because A' = 0. 




FIG. 3. Typical cuprate superconductor with Cu02 con- 
duction planes. 

We have so far identified several quantum group sym- 
metric models (with and without phonons) and have 
achieved a better understanding of ffext's superconduct- 
ing quantum symmetry. There are however still open 
questions: (i) Does a new model exist that is equivalent 
to Hq.sym in 1-D but can also be formulated on higher di- 
mensional lattices without breaking the symmetry? This 
would be important for realistic models, see Fig. 0. (ii) 
Are there models with symmetry away from half-filling? 
(iii) What is the precise relation between models with 
classical and quantum symmetry in this setting? 

As we shall see the answer to the last question also 
leads to the resolution of the first two. Without loss of 
generality (see argument given above) we will focus on 
one pair of next-neighbor sites in the following. We shall 
present two approaches that supplement each other: 

a. Generalized Lang-Firsov transformation We recall 
that the Hubbard model with phonons (with classical 
symmetry) can be transformed into the standard Hub- 
bard Hamiltonian in two steps: A Lang-Firsov transfor- 
mation changes the model to one with vanishing local 
phonon coupling and a mean field approximation removes 
the phonon operators from the model by averaging over 
Einstein oscillator eigenstates |2y]. There exists a sim- 
ilar transformation that relates the extended Hubbard 
model (with quantum symmetry) to the standard Hub- 
bard model: 



Hr 



H„ 



H, 



Hub- 



(We have already seen the first step of this transforma- 
tion above in (|lO|).) It is easy to see that the hopping 
terms of i/„_sym and i^Hub have different spectrum so 



the transformation that we are looking for cannot be an 
equivalence transformation. There exists however an in- 
vertible operator M, with MM* = 1 + (a^ - 1)^, C^ = C 
{i.e. similar to a partial isometry), that transforms the 
coproducts of the classical Chevalley generators into their 
Jinibo-Drinfel'd quantum counterparts 



MA,{H),M* = Ag{H),, 
and the standard Hubbard Hamiltonian into H„ 



MHii^^r 

This operator M is 



H„ 



M = 1 ® 1 + (a - l)e + /?/, 



(16) 



(17) 



(18) 



with / = X- ® X+ - X+ (E>X-,( = -f = i(ij2 (g> 
H^ - Hs® Hs) and a ± P = q^^ . With this knowledge 
the proof of the quantum symmetry of H^xt is greatly 
simplified. 

b. Quantum vs. Classical Groups — Twists A system- 
atic way to study the relation of quantum and clas- 
sical symmetries was given by Drinfel'd [gl|. He ar- 
gues that the classical U(g) and g-deformed Ug(g) uni- 
versal enveloping algebras are isomorphic as algebras. 
The relation of the Hopf algebra structures is slightly 
more involved: the undeformed universal enveloping al- 
gebra U(g) of a Lie algebra, interpreted as a quasi- 
associative Hopf algebra whose coassociator is an invari- 
ant 3-tensor, is twist-equivalent to the Hopf algebra Uq (g) 
(over [[Ing]]). 

All we need to know here is that classical (Ac) and 
quantum (Ag) coproducts are related via conjugation 
("twist") by the so-called universal J^ e Ug(su2)®^: 



Ag{x) ^ TA,ix)T-\ 



(19) 



(For notational simplicity we did not explicitly write the 
map that describes the algebra isomorphism of U(su2) 
and Uq(su2) but we should not be fooled by the appar- 
ent similarity between (|l^ and (|l^): The algebra iso- 
morphism does not map Chevalley generators to Jimbo- 
Drinfel'd generators and M is not a representation oi T.) 
The fundamental matrix representation of the univer- 
sal T for SU(N) is an orthogonal matrix |23] 



^(S2 j-j^-j ^ ^ eii (g) eu -f cos V3 ^ e^. 
+ sin if 2J (sij ® gji — Cji (E 

i<j 



(20) 



where cos (^ ± sin t^p = y^2q^''-/{q + q ^) , i,j — 1 . . .N 
and By are NxN matrices with lone "1" at position {i,j). 



The universal J-' in the superconducting spin-i represen- 
tation, i.e. essentially the N = 2 case with the Pauli 
matrices replaced by (H), is 



Fs = exp{ipf)s = ^ + cosip^ + sinip f 



(21) 



and ^-|-^ = ls(8)ls- We are interested in a representation 
of the universal J- on the 16-dimensional Hilbert space 
of states of two sites: 



J^=(em®Ps)®'(^)-exp((^/) 
= 1(K)1 - l=«ils + F,. 



(22) 



Note that the trivial magnetic representation em enters 
here even though we decided to study only deformations 
of the superconducting symmetry — Fs alone would have 
been identically zero on the hopping term and would 
hence have lead to a trivial model. 

We now face a puzzle: By construction F^^Hq_sy„iF 
should commute with the (global) generators of SUq(2)^ 
just like i^Hub- But F~^Hq_symF obviously has the same 
spectrum as i?g-sym so it cannot be equal to -ffjiub • There 
must be other models with the same symmetries. In fact 
we find a six-parameter family of classically symmetric 
models in any dimension. In the one-dimensional case 
twist-equivalent quantum symmetric models can be con- 
structed as deformations of each of these classical models. 
^Hub and Hq.sym are not a twist-equivalent pair but all 
models mentioned are related by generalized Lang-Firsov 
transformations. 

To close we would like to present the most general 
Hamiltonian with SU(2) x SU(2)/Z2 symmetry and sym- 
metric next-neighbor terms. (A group-theoretical deriva- 
tion and detailed description of this model is however be- 
yond the scope of this letter and will be given elsewhere.) 
The Hamiltonian is written with eight real parameters (/z, 
r, s, t, u, V, Re(z), Im(z)): 

Hsym = uy^n.in^x - f^'^n^^ + t y^ qLqj<t 
+ r ^ ni^rij-a + « X! '^i'^^J'^ 

Zfl U ^ — ^ 4. 4. 

+ — - — 2^ aiT"a"jT%l 
+ is-r) Y^ a\^a,^„a]_^aj„ 

+ v 22i ("-iT"'ii'^iT"'ji ^ rijf njjrijCT - nianj^rij^) 
+ ^ al_^aj^^ {z{ni„ - l)nj„ + z*ni„(nj^ - 1)) 
+ H.C. (23) 



For symmetry v = r + s + u~-2fi must hold. One pa- 
rameter can be absorbed into an overall multiplicative 
constant, so we have six free parameters. The first three 
terms comprise the standard Hubbard model but now 
without the restriction to half-filling. The filling factor 
is fixed by the coefhcient of the pair hopping term (6*'* 
term). The number e in the denominator of this coeffi- 
cient is the number of edges per site. For a single pair 
of sites e = 1, for a one-dimensional chain e — 2, for a 
honeycomb lattice e = 3, for a square lattice e = 4, for a 
triangular lattice e = 6 and for a _D-dimensional hyper- 
cube e = 2D. For a model on a general graph e will vary 
with the site. The 4*'' and 5*^ term describe density- 
density interaction for anti-parallel and parallel spins re- 
spectively. The balance of these two interactions is gov- 
erned by the coefhcient of the spin- wave term {7*^ term). 
The last term is a modified hopping term that is reminis- 
cent of the hopping term in the t-J model with hopping 
strength depending on the occupation of the sites; af- 
ter deformation this term is the origin of the non-trivial 
quantum symmetries of -ffsym ■ 

The known and many new quantum symmetric Hub- 
bard models can be derived from -ffgym by twisting as 
described above. While the deformation provides up to 
two extra parameters for the quantum symmetric mod- 
els the advantage of the corresponding classical models 
is that they are not restricted to one dimension. There 
are both classically and quantum symmetric models with 
symmetries away from half filling. 

The way the filling and the spin-spin interactions ap- 
pear as coefficients of the pair-hopping and spin-wave 
terms respectively looks quite promising for a physical 
interpretation. Due to its symmetries iJsym should share 
some of the nice analytical properties of the standard 
Hubbard Hamiltonian and could hence be of interest in 
its own right. 
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